A critical test of the mode-coupling theory of the glass transition 
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The mode-coupling theory ol the glass transition predicts the time evolution of the intermediate 
scattering functions in viscous liquids on the sole basis of the structural information encoded in two- 
point density correlations. We provide a critical test of this property and show that the theory fails 
to describe the qualitatively distinct dynamical behavior obtained in two model liquids characterized 
by very similar pair correlation functions. Because we use 'exact' static information provided by 
numerical simulations, our results are a direct proof that some important information about the 
dynamics of viscous liquids is not captured by pair correlations, and is thus not described by the 
mode-coupling theory, even in the temperature regime where the theory is usually applied. 
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I. INTRODUCTION 

Near the glass transition, the viscosity of a supercooled 
liquid is very sensitive to small temperature changes, 
and structural relaxation is likely a thermally activated 
process over an effective free-energy barrier, E, virhich 
grows upon cooling, r]{T) ^ r]aoexp{E{T)/kBT). From 
a theoretical perspective, one wishes to understand and 
provide accurate predictions for E{T) in terms of the 
microscopic interactions between particles in the liquid. 
However, even for the simple case of spherical particles 
with short-ranged, pairwise interactions, establishing a 
link between microscopic interactions and structural re- 
laxation remains an open theoretical challenge [l|. 

In this article, we assess one of the theoretical ap- 
proaches of the glass transition, the mode-coupling the- 
ory (MCT). The theory emerged in the mid-80s in the 
context of the theory of simple liquids and provides pre- 
dictions for the time evolution of density autocorrelation 
functions from the knowledge of the interaction between 
the particles [4]. Actually, the latter enters the theory 
only through two-point static density correlation func- 
tions. MCT therefore works under the hypothesis that 
knowledge of the pair correlation function g{r) of the 
fluid (or its Fourier transform S'(g), the static structure 
factor) is sufficient to predict the dynamical evolution. 
We focus here on this hypothesis and do not address the 
(known) failure of MCT to describe activated processes 
in deeply supercooled liquids. 

It is obvious that the evolution of the pair correlation 
function has some influence on the dynamics. Accord- 
ingly, some of the predictions that arc obtained within 
MCT compare reasonably well to numerical or experi- 
mental work. For instance, MCT provides guiding in- 
formation in those cases where the changes in the static 
pair correlation functions brought about by adding ei- 
ther very short-range attractions [3| , some asphericity Q 
or softness [5| to a dense system of hard spheres appear 
mostly responsible for the changes in dynamics reported 



for these systems. 

There are however also indications in the literature 
that there may be relevant structural information which 
is not captured by two-point functions. Broadly speak- 
ing, it is often stated that some form of 'local order' 
grows upon decreasing the temperature in supercooled 
liquids [1[. This order can be for instance (frustrated) 
icosahedral order [6| or the less well characterized 'amor- 
phous order' encoded in high-order 'point-to-set' correla- 
tion functions [3]. In both cases, g{r) is essentially blind 
to this growing static order, and the effect of the lat- 
ter on the dynamics (if any) would be completely missed 
by MCT. Similarly, it is by now understood that struc- 
tural relaxation near the glass transition occurs in a spa- 
tially correlated manner, with a dynamic lengthscale that 
grows upon cooling, with again no corresponding change 
in the pair correlation function [8| . 

Recently, we found that two model liquids with very 
similar structure at the level of pair correlation func- 
tions display distinct behavior upon cooling [9|. More 
specifically, we studied a binary Lcnnard- Jones (LJ) mix- 
ture [10| and the corresponding repulsive 'WCA mixture' 
in which attractive forces are truncated beyond the min- 
imum of the pair potential [ll[ , and we showed that the 
presence or absence of the attractive forces determines 
to a great extent the apparent energy barrier to struc- 
tural relaxation and its evolution with temperature and 
density. Therefore, these two model systems represent 
a benchmark on which the above ideas can be tested in 
much detail. 

In a previous work, Voigtmann used liquid-state the- 
ory to obtain the pair correlation function of the Lenard- 
Jones and the associated repulsive WCA models in the 
monoatomic case [12|- He then introduced the results 
in the MCT equations and found that the two liquids 
were essentially sharing the same dynamics, with nearly 
identical MCT critical temperatures. In our simulations 
on the other hand, we use a binary mixture, and we de- 
tect small differences between the structures of the two 



models. Thus, wc directly insert the 'exact' (numerically 
determined) structure factors into the MCT equations 
to test with a greater sensitivity the response predicted 
by MCT to small structural changes. Our main conclu- 
sion is that MCT fails to describe the difference in the 
dynamical behavior of the two liquids. More generally, 
this conclusion also suggests that microscopic theories of 
the glass transition based on pair correlation functions 
cannot accurately predict the evolution of the structural 
relaxation time. 

In Sec. |lT]we describe our models and method of anal- 
ysis. In Sec, mil we present our results for the dynamics. 
Section HVl concludes the paper. 



II. MODELS, METHODS AND STATIC 
RESULTS 



A. Static structure 



We consider the binary mixture of Lennard- Jones par- 
ticles introduced by Kob and Andersen [101 > ^^ '^^^^ ^^ i^^ 
WCA truncation [il|. The system is a 80:20 mixture of 
A:B particles interacting with the following interatomic 
pair potential between species a and /3, with a, /3 = A, B: 



'Jafi{r) = 4ea,3 
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In this expression, the cutoff 



distance r^o is equal to the position of the minimum of 
Va:j3 (r) for the WCA potential and to a conventional cut- 
off of 2. Sctq/j for the standard LJ model; Cap is a constant 
such that Vap (^qa) = 0. The value of the parameters aap 
and £„« were designed to obtain good glass-forming abil- 
ityM- 

We perform molecular dynamics simulations of both 
systems in the NVE ensemble after proper equilibration 
at a chosen state point characterized by a density p and 
a temperature T. We use N = 900 — 1300 particles (de- 
pending on the density) and work with periodic boundary 
conditions. A broad range of density has been consid- 
ered with p varying from 1.2 to 1.6. Lengths, tempera- 
tures and times are given in units of uaa, ^AA/ks, and 
(TOo-^^/48eAA)^^^ respectively. In line with the WCA 
theory [11|, the two liquid models are compared at the 
same (p, T) state points. Their pressures then differ, with 
the attractive interaction roughly providing a tempera- 
ture independent negative shift. 

From the numerical simulations we obtain both static 
and dynamic properties. We use static information as 
input for the MCT analysis described below, while the 
dynamical data are used to test the accuracy of the theo- 
retical results. For the dynamics, we record the evolution 
of the time dependence of the self part of the intermedi- 
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FIG. 1: Partial static structure factors SsBiq), SAAiq), and 
SAB{q) (from top to bottom at small g) at p = 1.2 near the 
theoretically determined Tc of MCT (top) and close to the 
Tc obtained from fitting the simulation data (bottom). The 
pair structures of the LJ and WCA liquids are very close, 
and the small deviations seem to increase somewhat when T 
decreases. 
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with qcTAA — 7.2, which corresponds to the position of 
the peak of the total static structure factor at p = 1.2, 
see Fig. [TJ In this expression, Na denotes the number of 
particles of species a, and rj{t) is the position of particle 
J belonging to species a at time t. From the decay of 
F^{q,t) for the majority species a = ^, we obtain the 
relaxation time t{p,T) which we conventionally define 
as Ff{q,T) = 1/e. Some of the data for t{p,T) were 
presented in Ref. Q. 

At the structural level, we measure the partial struc- 
ture factors Sapiq) which are needed as input for the 
MCT calculations. They are defined as 



Sap{q) = 
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^Nj^f 
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which can be written in a more compact form as a 2 x 2 
matrix S{q) whose matrix elements are the partial struc- 
ture factors Sapiq)- 




FIG. 2: Same as Fig. [T] for p = 1.6. Here the agreement 
between LJ and WCA mixtures is nearly perfect at all tem- 
peratures. 



We present a set of representative results for the par- 
tial structure factors in Fig. [l]for p = 1.2 and in Fig. [2] 
for p = 1.6. The former corresponds to the canonical 
density at which the Kob-Andersen mixture is usually 
studied [l3|- For each density, we present the data at 
two temperatures, corresponding to the location of the 
MCT transition determined either theoretically or by fit- 
ting the simulation data. Generically, the latter is lower 
than the former. 

These data confirm the well-known fact that the LJ 
and WCA models have a pair structure which is very sim- 
ilar, even in the case of a binary mixture. This supports 
the idea that the contribution of the attractive forces at 
the level of the two-point structure can be considered as 
a small perturbation and treated as such in theoretical 
liquid-state calculations [ll|. Looking more closely, we 
can nevertheless detect some small differences between 
the two systems. These differences increase slightly when 
one lowers the temperatures at constant density, as can 
be concluded for instance by comparing the two panels 
in Fig. [T] The agreement between the pair structures of 
the two liquids becomes excellent when density increases, 
see Fig. [2] This is physically reasonable since the steep 
repulsive core of the pair potential in Eq. ([T|), which is 
the same for both models, should play a predominant 
role under those thermodynamic conditions. 

We have obtained S{q) for a large number of state 
points at the densities p = 1.2, 1.4 and 1.6. We have then 



proceeded in two steps. First, we have directly inserted 
those data into the MCT dynamical equations presented 
below to obtain a rough determination of the location of 
the MCT transition. Then, we have used a smooth in- 
terpolation procedure to obtain the partial structure fac- 
tors in a small temperature window near the transition 
because it is not possible to resolve the structure directly 
from simulations with the accuracy needed to character- 
ize the MCT transition. This interpolation procedure is 
however quite innocuous, as the structure evolves in a 
mild and continuous manner in the temperature regime 
where the interpolation is performed. 



B. Mode-coupling theory 

The mode-coupling theory of the glass transition 
was originally derived to describe the dynamics of New- 
tonian systems [l3|. (A somewhat different derivation 
was also provided in the framework of the nonlinear fluc- 
tuating hydrodynamics [ij].) An analogous theory was 
later derived for Brownian systems |l5j . Here, we briefly 
present the latter version of MCT which, for convenience, 
we use in the following. 

The starting point of the theory is an exact equation 
for the time evolution of the intermediate scattering func- 
tions F(q; t) in terms of the so-called irreducible memory 
function. 



d_ 
dt 



Fiq,t) = -Doq^S-\q)Fiq,t) 



-J dt'M{q,t-t')-^Fiq;t'). 



(4) 



In this equation. Do is the diffusion coefficient of an iso- 
lated Brownian particle, and F(q, t) is a matrix whose ele- 
ments are the intermediate scattering functions F^^{q, t) 
defined as 



F'^P{q,t) 
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and is such that F{q,t ~ Q) — S{q). The matrix of mem- 
ory functions, M(q,i), can be expressed, in the mode- 
coupling approximation, in terms of the intermediate 
scattering functions F and structure factors S. The cor- 
responding (lengthy) expression can be found in Ref. 
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Equation (|4]) allows us to evaluate the time dependence 
of the intermediate scattering functions. The only input 
required is the static structure factor S{q). It is easy to 
see that the natural time unit for our system is a^ /Dq. 
In the following, all times are given in terms of this unit. 

Numerical resolution of Eq. (|3]) is difficult because one 
needs to describe the evolution of the intermediate scat- 
tering function on very widely separated time scales. The 
commonly used algorithm was first described in Ref. [1^1 ; 
here, we use the implementation described in great de- 
tail in Ref. [131. Briefly, the basic steps of the algo- 
rithm are as follows. The integro-differential equation is 



discrctizcd and solved numerically for 2Ns steps with a 
finite time step of dt. After 2Ng steps are completed, 
the time step is doubled and the results from the initial 
2Ns steps are mapped onto a new equally spaced set of 
Ns values for the quantities needed to continue the nu- 
merical algorithm. This mapping includes the integrals 
as well as the intermediate scattering functions. Then, 
the numerical algorithm is restarted with the new time 
step and proceeded for another Ng time steps, and the 
mapping is performed again. This procedure is followed 
until a convergence condition is satisfied. In the present 
work we used 200 equally spaced wave- vectors with spac- 
ing S = 0.2, the first wave-vector being at fcmin = S/2, 
and the largest wavevector at fcmax = 39.9. 



III. DYNAMICAL BEHAVIOR 

A. Intermediate scattering functions 

The solution of the MCT dynamical equations pro- 
vides the time evolution of the intermediate scattering 
functions F"^(g,i). We present a selection of data for 
the majority component, a ^ /3 = A, at density p = 1.2 
in Fig. [3l These data have the standard shape obtained 
within MCT, whereby the time decay streches over many 
decades of time when approaching the MCT singular- 
ity and occurs in two widely separated timescales. From 
now on, we focus on the slower of these, corresponding to 
the structural (alpha) relaxation of the system. Its time 
dependence is empirically well described by a stretched 
exponential form, 



Fiq,t) ^ cM-itM% 



(6) 



as commonly found in supercooled liquids. The data in 
Fig. [3] suggest that the stretching exponent /3 is essen- 
tially the same at all temperatures (and so MCT data 
obey a time-temperature superposition property) and is 
very close for both LJ and WCA hquids, with /3 « 0.87. 
We find similar values and agreement for all densities up 
to p = 1.6. For p = 1.2 and q = 7.25, Kob and Andersen 
report the value /3 = 0.85 flSj . 

The intermediate scattering functions of the two mod- 
els superimpose at very high temperatures, T = 3.0, 
when the structure factors also show very good agree- 
ment. However, when temperature is decreased, the dy- 
namics of the WCA model becomes faster than that of 
the LJ model, see the data for T = 1.0 in Fig. [3 The 
difference becomes more dramatic as the mode-coupling 
singularity is approached and long relaxation times are 
obtained. Although this is reminiscent of our numerical 
findings on the role of attractive forces in supercooled 
liquids Q , we shall see below that MCT only marginally 
accounts for this effect. 

As in numerical simulations, we determine the struc- 
tural relaxation time r(r, p) through the relaxation 
F°'^{q,T)/Saj3{q) = 1/e. This relaxation time depends 




FIG. 3: Normalized intermediate scattering function for the 
majority component for p = 1.2 and q = 7.1, for both LJ 
(full lines) and WCA (dashed lines) models obtained from 
the numerical solution of Eq. (Q. The temperature decreases 
from left to right. LJ: T = 3.0, 1.0, 0.908, 0.899, 0.8975, 
0.8972, 0.8971. WCA: T = 3.0, 1.0, 0.77, 0.746, 0.7425, 
0.7421, 0.74198. Note that the two sets of temperatures are 
not identical. The lower temperatures in each set provide a 
good estimate of the Tc values for each model. 



on the wavevector q, the chosen species a and f3. It is 
a priori different from the time extracted from simula- 
tions where only the incoherent part of the intermediate 
scattering function is considered. However, a simplifying 
feature of MCT calculations is that the evolution of r is 
the same for all observables near the transition. Thus, 
it does make sense to focus on a given wavevector (we 
choose the maximum of the first diffraction peak) , to fo- 
cus on the majority component of the liquid (a = /? = A), 
and to compare these results to the simulation data for 
which the incoherent scattering function is measured. 

The only difference between the two sets of data shown 
in Fig. [3] lies in the temperature regimes they cover 
since the MCT transition of the WCA model occurs at a 
slightly lower temperature than that of the full L J model 
including the attractive component of the potential. We 
now turn to a more precise discussion of this difference. 



B. Relaxation times 

We now enter the core of our analysis and study in 
detail the evolution with T and p of the structural re- 
laxation time, as predicted by MCT, and compare these 
theoretical results to those from computer simulations. 

We first present in Fig. 3] our results in an Arrhe- 
nius plot, showing the evolution of t{T, p) on a log-scale 
against the inverse temperature 1/T, for both the LJ and 
the WCA models, as measured in MCT calculations and 
in simulations. We show results for p = 1.2 and p = 1.6. 

In this Arrhenius representation, the results of MCT 
calculations appear almost 'vertical', because MCT pre- 
dicts the existence of a critical singularity Tc at which 
the relaxation time diverges algebraically |2|. This di- 



vergcncc is not seen in computer simulations and so the 
simulation data look qualitatively different. 

A second striking observation from Fig.|4]is that MCT 
clearly overestimates the temperature regime where slow 
dynamics sets in by about 100 %. Thus, there is no range 
of temperature where the theoretical calculations follow 
the simulation results, even if a vertical shift (correspond- 
ing to a trivial rcscaling of the microscopic timcscale) is 
allowed. 

At high temperatures, a good empirical representation 
of both MCT and simulation results is obtained by in- 
troducing an Arrhenius-like temperature dependence. 



>exp 



Eoc 



(7) 



Although there is no specific theoretical foundation for 
such a fit, it often accounts quite well for high tempera- 
ture results, even in experiments [1^. When fitting the 
data to Eq. (7), we obtain the values Eoo ~ 5.5 and 
4.2 for the LJ and WCA systems within MCT, while 
the simulation data are characterized by Eoo ~ 2.55 and 
2.0. The different values found for the two models in the 
high-temperature regime reflect the fact that attractive 
forces already play an important quantitative role in this 
regime, despite the very little changes observed in static 
quantities. This effect was reported before in the case of 
simple hquids [20|-[2^. 

The role of attractive forces at high temperatures is 
partially captured by our MCT calculations, since we find 
theoretically a 25 % change in Eoo between the two mod- 
els. However, given that the predicted values are off by a 
factor of 2, this agreement could well be fortuitous. At a 
density p = 1.6, we find E^o f* 12 from simulations of the 
two systems, while MCT again overestimates this value 
and predicts E^o ~ 22.5 for both models. So, at high 
temperatures and large densities, the difference between 
LJ and WCA mixtures vanishes both in simulations and 
MCT calculations, but the activation energy predicted 
by the latter is still in error by a factor close to 2. It is 
interesting to note that there is nearly the same factor 2 
in error for Eoo and for Tc (see below), which suggests 
that MCT does not merely 'break down' at low T but is 
in fact always quite inaccurate in its prediction for the 
temperature evolution of the structural relaxation time, 
even in the high-temperature liquid. 

The strong disagreement between MCT results and 
simulations is well-known Q and is two-fold: (i) MCT 
overstimates the temperature regime where slow dynam- 
ics takes place; (ii) the predicted algebraic divergence of 
the relaxation time is not really observed. When fitting 
data to MCT predictions, one can thus only find a mod- 
est temperature regime where an algebraic divergence de- 
scribes the data and one must simultaneously resort to 
some sort of rescaling or shift of the control parameters. 
Since this gives considerable freedom in the data treat- 
ment, we adopt the following procedure, shown in Fig.U) 
At a given density, we first concentrate on the LJ model. 
We then rescale r and T by adjustable time and tem- 
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FIG. 4: Comparison of the relaxation times obtained by MCT 
(open symbols) and simulations (filled symbols) at p = 1.2 
(upper panel) and p = 1.6 (lower panel). We both show the 
bare MCT data and the data rescaled to provide the best col- 
lapse of the L J data. The rescaling factor for the temperature 
is 2 for p = 1.2 and 1.9 for p — 1.6. This same factor is then 
applied with no further changes to the WCA model. 



perature scales in order to obtain the 'best' collapse of 
theory and simulations. We obtain a temperature scaling 
factor of 2 for p = 1.2, and 1.9 for p = 1.6, which mirrors 
the factor 2 in Eoo found above at high temperatures. 
Once this rescaling is performed for the LJ model, we 
use the same scaling factor for temperature in the WCA 
model. The results are presented in Fig. |4l While the 
rescaling works well if T is not too low in the LJ system, 
it performs very poorly for WCA, even when density is as 
large as p = 1.6. This implies that the different dynam- 
ical behaviors observed in simulations for LJ and WCA 
models [9| is only very qualitatively reproduced by MCT 
and is considerably underestimated by the theory at a 
quantitative level. 

This last statement can be made more precise by de- 
termining the location of the MCT singularity. We first 
determine Tc within the theory by following the growth 
of r near the transition, which is well described by an 
algebraic divergence: 

T^iT-Tc)-^. (8) 

The value of the critical exponent 7 depends in principle 
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FIG. 5: MCT plot of the relaxation times in which the alge- 
braic divergence in Eq. ([8]) appears as a straight line. Top: 
theoretical curves fitted with the same 7 = 2.4 for p — 1.2, 
1.4, and 1.6 (from bottom to top). Bottom: same for nu- 
merical simulations, but note the difference in scales and the 
deviations occurring at low temperatures. Data at different 
densities are vertically shifted, for clarity. 





1.2 


1.4 


1.6 


LJ - MCT 


0.8971 


1.8677 


3.528 


LJ - Simulation 


0.435 


0.93 


1.76 


WCA - MCT 


0.7419 


1.7707 


3.489 


WCA - Simulation 


0.28 


0.81 


1.69 



TABLE I: Theoretical critical temperatures Tc for the LJ and 
WCA models at different densities, and comparison with es- 
timates from fits to the simulation data. 



should appear as a straight line. The fit holds over about 
2-3 decades in relaxation time, as commonly found in 
the analysis of numerical and experimental data [3 , and 
deviations at low temperatures occur similarly for both 
models at all densities. As noted before [131, it should 
be remarked that the algebraic fit to simulation data is 
obtained in a range of (T — Tc)/Tc where the theoreti- 
cal predictions have not yet entered the critical regime, 
which casts some doubts on the consistency of the fitting 
procedure. 

The values of the critical temperatures obtained by 
fitting the simulation data are also reported in Table 
m As found theoretically, the difference between the 
LJ and WCA models is smaller at larger density, when 
the structure factors become more similar (see Fig. [2|). 
However, the difference in critical temperatures decreases 
from 36 % at p = 1.2 to 4 % at p = 1.6 and is much larger 
than the MCT theoretical predictions at all densities (re- 
spectively 17 % and 1 %). Thus, the difference in the 
dynamical behavior of the two models is only marginally 
captured by MCT calculations and is quantitatively con- 
siderably underestimated. 



on the studied model. In Fig. [SJ we show however that 
the value 7 = 2.4 actually describes both the LJ and the 
WCA models for densities p = 1.2, 1.4 and 1.6. The 
values of the critical temperatures are reported in Table 
|T) Tc increases by a factor « 4 when density increases 
from 1.2 to 1.6. As density increases, the difference in 
critical temperatures between the LJ and WCA models 
becomes smaller, decreasing from 17 % at p = 1.2 to 1 % 
at p = 1.6 (within the theory). 

While the determination of Tc is unambiguous in MCT 
calculations, this is not the case in simulations where the 
MCT power law cannot be followed for arbitrarily large 
relaxation times. In practice, this means that the values 
of Tc and 7 obtained in simulations are strongly corre- 
lated and therefore depend on the chosen temperature 
range for fitting. Thus, we decided to constrain our data 
analysis by imposing the value of the exponent 7 obtained 
in the theoretical calculations. We then determined the 
value of Tc that fits the data best with this exponent. The 
quality of the fits can be judged in Fig. [51 where the data 
for t{T, p) are plotted as a function of the reduced vari- 
able {T — Tc)/Tc in a log-log representation where Eq. ([5|) 



IV. DISCUSSION 

The structure of simple liquids is usually described in 
terms of pair correlation functions, and the development 
of analytical theories to predict their evolution with den- 
sity or temperature from the knowledge of the interaction 
between particles has been a major theoretical achieve- 
ment [23 ■ Whether knowledge of ^(r) is enough to char- 
acterize the structure and predict the dynamics of viscous 
liquids has however been an open question in the field 
of the glass transition. This is a central issue for mode- 
coupling theory which makes dynamical predictions from 
the sole knowledge of two-point density correlations. 

To address this issue we have used the pair correla- 
tion functions taken from molecular dynamics simula- 
tions to obtain MCT predictions for the relaxation time 
of two model liquids characterized by similar pair struc- 
tures, but very distinct dynamics [a]- We have found that 
MCT is essentially blind to these dynamical changes, be- 
ing unable to account for them in a quantitative way. In 
this case at least, the necessary structural input for the 
dynamics of the viscous liquids is not merely encoded 
in two-point density correlation functions. Accordingly, 



MCT is bound to yield quantitatively inaecurate pre- 
dictions. As mentioned in the introduction, there still 
remains room for applying MCT, as there exists physi- 
cally relevant glassy phenomena associated to important 
changes in the structure that are captured by two-body 
correlations. 

In a similar vein, the present work has been an op- 
portunity to test the idea, popularized in particular in 
the context of MCT, that even small changes in the pair 
structure may have dramatic consequences in the dynam- 
ics. For the case under study, we have found that the 
small differences seen between the structure factors of 
the two liquid models only produce, within MCT, mi- 
nor dynamical differences, thus confirming Voigtmann's 
results [12|- As already stated, these MCT predictions 
dramatically underestimate the changes observed in the 
simulations. An additional factor explaining this dis- 
crepancy is the fact that MCT overestimates by a large 
amount the location of the putative singularity (empiri- 
cally estimated by fitting simulation data). As a result, 
MCT predictions are based on structure factors measured 
at relatively high temperatures where differences between 
the WCA and LJ models arc even less pronounced: com- 
pare both panels in Fig. [T] 

Finally, our work does not address whether the failure 
exposed here of the 'g{r) determines the dynamics' MCT 
motto that is shown here is connected to the failure of 
MCT to describe activated dynamics at low temperature. 
For instance, in their attempt to incorporate activated 
processes in the MCT framework, Schwcizer and cowork- 
ers [2J] manage to avoid the MCT algebraic divergence 



of Eq. ([8]). However, in their approach, it is still the pair 
correlation function which determines the amplitude of 
the free energy barriers that have to be crossed dynam- 
ically. In this case, one also expects the theory to miss 
the difference of dynamical behavior between the L J and 
WCA liquid models, despite the presence of activated 
processes. 

To conclude, our work shows that structural informa- 
tion not incorporated in two-body density correlations 
likely play an important quantitative part in driving the 
slowdown of dynamics in viscous liquids. This raises se- 
rious doubts on the ability of a number of analytical ap- 
proaches to make quantitative predictions in supercooled 
liquids. In this paper, we have only focused on MCT 
and we leave for a future publication a more general dis- 
cussion of the consequences of this finding. Our results 
also motivate further research to detect more complicated 
forms of 'hidden', 'amorphous', or 'frustrated' order in 
liquids approaching the glass transition. 



Acknowledgments 

We wish to thank E. Flenner and G. Szamel for gener- 
ously making their numerical code to solve MCT equa- 
tions for a binary mixture [17[ available to us, thereby 
making this work possible. We also thank W. Kob for 
useful advice. We acknowledge partial support from the 
ANR Dynhet. 



[1] P. G. Debenedetti and F. H. Stillinger, Nature 410, 259 [12 

(2001). [13 

[2] W. Gotze, Complex dynamics of glass-forming liquids: A 

mode-coupling theory (Oxford University Press, Oxford, [14 

2008). 
[3] K. N. Pham, A. M. Puertas, J. Bergenholtz, S. U. Egel- [15 

haaf, A. Moussaid, P. N. Pusey, A. B. Schofield, M. E. [16^ 

Gates, M. Fuchs, and W. C. K. Poon, Science 296, 5565 

(2002). [17 

[4] M. Letz, R. Schilling, and A. Latz, Phys. Rev. E 62, 5173 

(2000). [18 

[5] L. Berthier, E. Flenner, H. Jacquin, and G. Szamel, Phys. 

Rev. E 81, 031505 (2010). [19 

[6] G. Tarjus, S. A. Kivelson, Z. Nussinov, and P. Viot, J. 

Phys.: Gondens. Matter 17, R1143 (2005). [20 

[7] A. Montanari and G. Semerjian, J. Stat. Phys. 125, 22 

(2006). [21 

[8] M. D. Ediger, Annu. Rev. Phys. Ghem. 51, 99 (2000). 
[9] L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 [22 

(2009). 
[10] W. Kob and H. G. Andersen, Phys. Rev. Lett. 73, 1376 [23 

(1994). 
[11] J. D. Weeks, D. Ghandler, and H. G. Andersen, J. Ghem. [24 

Phys. 54, 5237 (1971). 



T. Voigtmann, Phys. Rev. Lett. 101, 095701 (2008). 
U. Bengtzelius, W. Gotze, and A. Sjolander, J. Phys. G: 
Sohd State Phys. 17, 5915 (1984). 

S. P. Das, G. F. Mazenko, S. Ramaswamy, and J. J. Toner 
Phys. Rev. Lett. 54, 118 (1985). 

G. Szamel and H. Lowen, Phys. Rev. A 44, 8215 (1991). 
M. Fuchs, W. Gotze, I. Hofacker, and A. Latz, J. Phys.: 
Gondens. Matter 3, 5047 (1991). 

E. Flenner and G. Szamel, Phys. Rev. E 72, 031508 
(2005). 

W. Kob and H. G. Andersen, Phys. Rev. E 52, 4134 
(1995). 

G. Alba-Simionesco, D. Kivelson, and G. Tarjus, J. 
Ghem. Phys. 116, 5033 (2002). 

J. Kushick and B. J. Berne, J. Ghem. Phys. 59, 3732 
(1973). 

D. P. Dean and J. N. Kushick, J. Ghem. Phys. 76, 621 
(1982). 

S. D. Bembenek and G. Szamel, J. Phys. Ghem. B 104, 
10647 (2000). 

J. P. Hansen and I. R. McDonald, Theory of Simple Liq- 
uids (Elsevier, Amsterdam, 1986). 

See K. S. Schweizer, J. Ghem. Phys. 127, 164506 (2007) 
and refs. therein. 



